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Abstract 

We report the first observation of color-suppressed B^ — > D^tt^ and D^*'^lo 
. decays and evidence of B^ — > 0*^7^^ and D^*'^r]. The branching fractions 

(X) '. are found to be B{W -^ D^tt^) = (2.9 toj ± 0-6) x 10-^ B{W -^ D^lu) = 

^ ; (1.7 t[]ito:l) X 10-^ and B{W -^ D*^uj) = (3.4 t};?±0.8) x lO"! The anal- 

ysis is based on a data sample of 21.3 fb^^ collected at the T{4S) resonance 



^D ■ by the Belle detector at the KEKB e^e collider. 
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The decay modes B^ -^ D^*^^X^, where X'^ is a hght neutral meson, proceed via an 
internal spectator diagram and are expected to be suppressed relative to the external dia- 
gram, since the color of the u antiquark produced by the weak current must complement the 
color of the c quark as shown in Fig. |T[ Studies of such color-suppressed decay modes can 
be used to test models of hadronic B meson decays and provide information on final-state 
interactions. Results for color-suppressed 5° -^ /)(*)0x'^ decays have been published by the 
CLEO collaboration [0]; however, only upper limits were obtained. 

In this paper, we report on a search for the color-suppressed B^ -^ D^X^, and D*^X^ 
decay processes, where the neutral meson X^ is either a 7r°, t], or uj. Charge conjugate modes 
are implicitly included in this paper. The data sample used in this analysis was collected 
with the Belle detector at KEKB . KEKB is a double storage ring with 8 GeV electrons 
and 3.5 GeV positrons colliding at a 22 mrad crossing angle. The data sample corresponds 
to an integrated luminosity of 21.3 fb~^ at the T(4S') resonance and contains 22.8 million 
BB pairs. 

Belle is a general-purpose detector containing a 1.5 T superconducting solenoid mag- 
net. Charged particle tracking covering 92% of the total center-of-mass (CM) solid angle is 
provided by a Silicon Vertex Detector (SVD) consisting of three concentric layers of double 
sided silicon strip detectors, and a 50-layer Central Drift Chamber (CDC). Particle identifi- 
cation is accomplished by combining the responses from an array of Silica Aerogel Cerenkov 
Counters (ACC) and a Time of Flight Counter system (TOF) with dE/dx measurements in 
the CDC. The combined response of the three systems provides at least 2.5(7 K/tt separation 
for laboratory momentum up to 3.5 GeV/c. Photons and electrons are detected in an array 
of 8736 CsI(T£) crystals (ECL) located inside the magnetic field and covering the entire 
solid angle of the charged particle tracking system. The 1.5 T magnetic field is returned 
via an iron yoke that is instrumented to detect muons and K^ mesons (KLM). The KLM 
consists of alternating layers of resistive plate chambers and 4.7 cm thick steel plates. 

For the light neutral meson X°, we use the vr*^ -^ 77, r] -^ 77, rj -^ tt^tt^tt^ and the u -^ 
7r+7r~7r° decay channels. Charged tracks are required to have impact parameters that are 
within ib5 cm of the interaction point along the positron beam axis and 1 cm in the transverse 
plane. We reject tracks that are consistent with electrons or muons. The remaining tracks 
are identified as pions or kaons according to a kaon to pion likelihood ratio. Candidate 7r° 
mesons are reconstructed from pairs of photons in the ECL that have an invariant mass 
within zbl6 MeV of the nominal tt" mass Q]. The vr'' daughter photons are required to have 
energies greater than 50 MeV. Both photons from the 77 — >■ 77 mode are required to have 
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E^ > 100 MeV and the energy asymmetry of the daughter photons, ^"^^_)_^^ , is required to 
be less than 0.8. We remove t] candidates if either of the daughter photons can be combined 
with any other photon with E^ > 100 MeV to form a tt^ candidate. Candidate t] mesons 
for the 77 mode are required to have an invariant mass within ±10.6 MeV/c^ (±2.5(t) of 
the nominal 7] mass; the corresponding requirement for the 7r"'"7r~7r° mode is ±3.4 MeV/c^. 
Candidate r] mesons are constrained to the nominal rj mass; the fit also constrains the tt+tt" 
pair from the Sir channel to a common vertex point. Candidate u mesons are tt+tt^tt^ 
combinations with an invariant mass within ±30 MeV/c^ of the nominal u mass value. The 
CM momentum of the 7r° from the u decay is required to be greater than 350 MeV/c to 
reduce the large combinatorial background from low energy photons. The vr+vr" pair from 
the uj decay is required to form a common vertex point within the beam interaction region 



taking into account the lifetime of the B meson. Detailed studies of tracking, tt" detection, 
and particle identification yield systematic errors in the detection efficiencies of 7.3% for 
prompt 7r°, 6.6% for rj, and 7.1% for uj mesons. 

For candidate D^ mesons we use the D° — ^ K^tt^, K^tt^tt^, and K^ti^ti^ti^ decay 
modes. The CM momentum of the 7r° from D'^ -^ K^ti^'k'^ decay is required to be greater 
than 300 MeV/c. The invariant mass of the D^ candidates are required to be within ±2.5a 
of the measured D^ mass, where a is the D^ mass resolution that varies between 5.5 and 13 
MeV/c^ depending on the decay mode. A mass and vertex constrained kinematic fit is then 
performed to the D^ candidates. Good D^ candidates are required to have an acceptable x^ 
value from the fit. D*^ candidates are reconstructed in the D*^ —>■ D^tt^ decay mode. For 
these 7r° mesons, the photon energy cut is reduced to 20 MeV. For D*° candidates, the mass 
difference. Am = M(D°7r°) -M(D°), is required to be within 2.5a (a = 0.8 MeV/c^) of the 
nominal mass difference. The systematic error on the D° meson reconstruction efficiency 
is studied using the B~ — >■ D^vr" data sample. The error is determined to be 15% from a 
comparison of the observed number of signal events relative to the expected yield assuming 
the PDG B^ -^ D^Ti^ branching fraction [Q. The systematic uncertainty in the detection 
efficiency for low momentum 7r*^s is found to be 10.7%. 

We combine D°s or D*^s with X° meson candidates to form 5° candidates. Two kine- 
matic variables are used to identify signal candidates, the beam-constrained mass M^c = 



(^beam)^ - (pg^)^ and the energy difference AE = E^^ - E^^^, where E'^^ and pf^ are 
the CM energy and momentum of the B^ candidate, and -E'^e'am ~ V^/2 — 5.290 GeV. The 
typical Mbc resolution is 3 MeV/c^; the AE resolution ranges from 17 to 25 MeV, depending 
on the decay mode. When more than one B'^ candidate is found in an event, the candidate 
with the minimum x^ is chosen, where x^ = Xdo + Xxo(+XAm)- Here Xdo is the x^ of the 
kinematic fit to the D^, x'xo is the x^ of the kinematic fit to either the 7r° or the r/. For 
the u, Xxo = {A{M^)/a{Mi^)y, where a is the measured resolution. For the D*°X° modes, 
XAmy defined as (A(Am)/a(Am))^, is included in the best candidate selection. 

The background from continuum e~^e~ -^ qq production is suppressed in the following 
ways. For the D^*^^t] final state, we apply cuts on the ratio of the second to zero-th Fox- 
Wolfram moments 0, i?2, and the angle between the thrust axis of the B candidate and 
the thrust axis of the rest of the event (cos^t)- For the D^*^^7i^ and D^*'>^u final states, 
we use the B fiight direction and a Fisher discriminant ^ containing several variables that 
quantify event topology . We also use the D*^ helicity angle for the D*^7t^ mode and the uj 
helicity angle and the cu decay amplitude for the D^*^'^uj mode |Q. Each of these variables is 
parameterized to form signal (S) and background (BG) probability density functions (PDF). 
The PDFs are multiplied to form a single likelihood £s(bg)) and then a cut is apphed on the 
likelihood ratio C^/i'^s + -^bg) to suppress the qq background. Signal PDFs are determined 
using Monte Carlo (MC) and background PDFs are obtained from Mf,c sideband data. The 
cut efficiencies are typically 70% and remove more than 90% of the qq background. The 
systematic error in the efficiency for this cut is determined by applying the same procedure 
to the 5" -^ D*^-K~ data sample. By comparing the cut efficiency between the data sample 
and MC, the systematic error is determined to be 5%. 

In addition to the qq background, we observe large background contributions from color 
favored B -^ D^*\mi)~ decays and cross-talk from D*^X^ to D°X° modes. The D*~^p~ 
mode has the same final state as D^u and D^rj. However, when the D*^ daughter pion 



is combined with the p~, the invariant mass rarely falls within the uo oi t] mass window. 
The D^*^^p~ final state contaminates the D^*^^Tr'^ mode if the p~ decays to a fast vr*^. This 
mode also contaminates the D^*^^t] channel if a photon from the fast tt" is combined with 
another photon to form an rj candidate. About half of these events are removed by explicitly 
reconstructing the D^*'^p^ final state. The contributions of these backgrounds in the ?y 
channel, as well as the feed-across from the Z)'^*^°7r° mode, is also minimized by the 7r° 
veto discussed above. We also check for background contributions from B -^ D^*^^ p'~ 



{p'~ — *■ cuTT") decays which have recently been observed by CLEO [|T0|. This two-body decay 
produces high momentum D^*^^s and cos that can fake signal events. Monte Carlo studies 
indicate that the remaining background events are shifted in AE by approximately the mass 
of the missing slow pion and thus can be distinguished from signal events by fitting the AE 
distribution. 

The AE distributions for the various D^*''^X^ decays are shown in Fig. ^ after applying all 
selection cuts and requiring Mbc to be between 5.272 GeV/c^ and 5.288 GeV/c^. The signal 



is modeled with a Crystal-Ball function |[TT]| with parameters obtained from MC. The back- 
ground functions include a combinatorial component and a color- favored component. The 
D^X^ modes also include a component for cross-talk from color- suppressed D*^X^ modes. 
The combinatorial component is taken to be a second order polynomial with parameters 
determined by the AE shape in the M^^ sideband (5.20 GeV/c^ < M^c < 5.26 GeV/c^). 
The shapes of the color-favored and cross-talk components are modeled by MC histograms. 
The signal and background normalizations are free parameters in each fit. 

Table H lists the signal yield, statistical significance, efficiency, and the branching frac- 
tion for each D^*^^X'^ mode. The systematic errors due to fitting are obtained by varying 
the parameters of the fitting functions within la of their nominal values. The change in 
the signal yield from each variation is added in quadrature to obtain the fitting systematic 



errors. These are typically 10%. The statistical significance is defined as w— 21n(£(0)/£inax) 
where £max is the likelihood at the nominal signal yield and C{0) is the likelihood with 
the signal yield fixed to zero. We observe signals for B^ —>■ D'^tt^, D^uo, and D*^uj decays 
with more than 4a significance. We find evidence of signals for 5° — *■ D*^7t^, D^rj, and 
D*^ri with more than 3a significance. For decay modes with significance less than 4, we 
give 90% confidence level upper limits (UL) on the signal yields [N^^) from the relation 

/o " C{Ns) dNs I J^ C{Ns) dNs = 0.9, where C{Ns) denotes the maximum likelihood with 
the signal yield fixed at A^s- The efficiencies for each decay mode are calibrated with con- 
trol data samples. The final systematic errors include the errors in fitting, reconstruction 
efficiency, background suppression cut efficiency, and the number of BB pairs. Assuming 
the number of B^B^ and B^B~ pairs are equal, we calculate the branching fractions for 
various decay modes given in Table |. The branching fraction upper limits are calculated 
by increasing A"g^^ and reducing the efficiency by their systematic errors. 

In summary, using 22.8 million BB events collected with the Belle detector, we report the 
first observations of color-suppressed B^ — > D^tc^ and D^*^^uj decays. We also find evidence 
for B^ — > D*^7r^ and D^*^^r] signals. All the color-suppressed modes have similar branching 
fractions with central values between 1.4 and 3.4 xlO~^, as shown in Table |. In general, 
the branching fractions are consistently higher than recent theory predictions [Q based on 
the factorization hypothesis. This may be accounted for by additional corrections to the 
factorization models, or by non-factorizable effects such as final state interactions. 
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TABLE I. The obtained signal yield, statistical significance, efficiency including the sub-decay 

branching fractions, branching fraction (B), and 90% confidence level upper limit (UL) for each 
;g0 ^ ^(*)Oj^o jgpj^y ^Q^g_ 
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FIG. 1. The external(left) and internal(right) spectator diagrams for B -^ Dtt decays. 
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FIG. 2. The AS distribution for (a) D^tt^, (b) D*^tt^, (c) D^t], (d) D*^r], (e) D^lo, and (f) D*^u;. 
The sohd hne shows the fitting result. The dashed hne shows the sum of the signal component 
and the combinatorial background component. The combinatorial component is shown separately 
as the cross-hatched area. 
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